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Stochastic  Lagrangian  models  for  the  velocity  following  a  fluid  particle  are  used  both  in  studies  of 
turbulent  dispersion  and  in  probability  density  function  (PDF)  modeling  of  turbulent  flows.  A 
general  linear  model  is  examined  for  the  important  case  of  homogeneous  turbulent  shear  flow,  for 
which  there  are  recent  direct  numerical  simulation  (DNS)  data  on  Lagrangian  statistics.  The  model 
is  defined  by  a  drift  coefficient  tensor  and  a  diffusion  tensor,  and  it  is  shown  that  these  are  uniquely 
determined  by  the  normalized  Reynolds-stress  and  timescale  tensors  determined  from  DNS.  With 
the  coefficients  thus  determined,  the  model  yields  autocorrelation  functions  in  good  agreement  with 
the  DNS  data.  It  is  found  that  the  diffusion  tensor  is  significantly  anisotropic — contrary  to  the 
Kolmogorov  hypotheses  and  conventional  modeling — which  may  be  a  low-Reynolds-number  effect. 
The  performance  of  two  PDF  models  is  also  compared  to  the  DNS  data.  These  are  the  simplified 
Lagrangian  model  and  the  Lagrangian  isotropization  of  production  model.  There  are  significant 
differences  between  the  autocorrelation  functions  generated  by  these  models  and  the  DNS  data. 
©  2002  American  Institute  of  Physics.  [DOI:  10.1063/1.1465421] 


I.  INTRODUCTION 

Homogeneous  turbulent  shear  flow  is  of  fundamental 
importance  in  the  development  of  models  for  inhomoge¬ 
neous  turbulent  flows.  Both  experiments1  and  direct  numeri¬ 
cal  simulations  (DNS)2  of  homogeneous  shear  flow  have 
been  performed  in  which  Eulerian  statistics  of  the  turbulence 
have  been  measured.  More  recently,  a  series  of  DNS  studies 
has  been  performed3 -5  in  which  Lagrangian  statistics  have 
been  obtained  by  tracking  a  large  number  of  fluid  particles. 
These  studies  clearly  have  direct  relevance  to  stochastic  La¬ 
grangian  models6  of  turbulence,  which  model  the  motion  of 
fluid  particles  as  diffusion  processes  (i.e.,  continuous  Markov 
processes).7  The  purpose  of  this  paper  is  to  show  the  connec¬ 
tion  between  the  Lagrangian  velocity  autocovariance  tensor 
obtained  from  DNS  and  stochastic  Lagrangian  models  for 
fluid  particle  velocity. 

Stochastic  Lagrangian  models  for  the  velocity  of  a  fluid 
particle  arise  in  two  different  contexts:  turbulent 
dispersion;8-10  and  probability  density  function  (PDF) 
models.1112,7  In  both  cases  the  general  form  of  the  models 
considered  (when  applied  to  homogeneous  turbulence)  can 
be  written  as  the  linear  stochastic  differential  equation  (SDE) 

duj= —AjjUj  dt  + B jjdWj ,  (1) 

where  du(t)  =  u(t  +  dt)  —  u(f)  is  the  infinitesimal  increment 
of  the  fluctuating  component  of  velocity  u(f)  following  the 
fluid  particle;  we  refer  to  A(  t )  as  the  drift  tensor;  B(  t )  is  the 
diffusion  coefficient;  and  dW(t)  is  the  infinitesimal  incre¬ 
ment  of  a  vector-valued  Wiener  process  which  has  the  prop¬ 
erties  (d W)  =  0,  {dW i  dWj)  =  dtSjj  .  Different  models  corre¬ 
spond  to  different  specifications  of  the  drift  tensor  A(f)  and 
diffusion  coefficient  B(f). 


For  statistically  stationary,  homogeneous  isotropic  turbu¬ 
lence  (with  no  mean  velocity  gradients)  the  only  sensible 
choice  of  coefficients  is 


and 


BU= 


2  m  ' 2  \ 1/2 


(2) 

(3) 


where  TL  is  the  Lagrangian  integral  timescale  and  u '  is  the 
turbulence  intensity  (i.e.,  the  rms  velocity  fluctuation).  Then, 
Eq.  (1)  reduces  to  an  independent  Langevin  equation  for 
each  component  of  velocity 

dt  1 2u'2\m 

dn,=  -Mi— +  |— )  dWi.  (4) 


This  model  dates  back  to  Taylor’s  1921  original  paper  on 
turbulent  dispersion.8  The  autocorrelation  function  given  by 
Eq.  (4)  is 

p(s)  =  (ui(t)uft  +  s))/M,2  =  exp(  —  \s\/Tl),  (5) 

which  agrees  well  with  DNS  data13  (except  at  small  values  of 

M  itl). 

The  central  issue  addressed  here  is  the  appropriate  speci¬ 
fication  of  A  and  B  in  homogeneous  turbulent  shear  flow. 
This  has  been  considered  in  the  context  of  turbulent  disper¬ 
sion  by  Sawford  and  Yeung.4'5  These  authors  compared  La¬ 
grangian  autocorrelations  predicted  by  two  dispersion  mod¬ 
els  to  DNS  data.  Both  of  these  models  take  B  to  be  isotropic. 

We  show  here  that  appropriate  values  of  A  and  B  can  be 
deduced  from  the  measured  Lagrangian  velocity  autocovari¬ 
ance,  and  that  the  resulting  model  is  in  good  agreement  with 
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the  DNS  data.  This  agreement  supports  the  nontrivial  con¬ 
clusion  that  the  Lagrangian  velocity  is  well  represented  by  a 
linear  diffusion  process  (except  over  small  time  intervals). 
The  deduced  value  of  B  is  significantly  anisotropic. 

The  performance  of  two  models  used  in  PDF  methods  is 
compared  to  the  DNS  data.  These  are  the  simplified  Lange- 
vin  model  (SLM)  and  the  Lagrangian  isotropization  of  pro¬ 
duction  (LIPM)  model.14 


II.  HOMOGENEOUS  TURBULENT  SHEAR  FLOW 


In  homogeneous  turbulent  shear  flow,  the  imposed  mean 
velocity  gradient  is 


3{U,) 


dXj 


—  SSj  [  Sj2 , 


(6) 


where  S  is  the  (constant)  imposed  mean  shear  rate.  The  tur¬ 
bulence  is  characterized  by  the  Reynolds  stress  tensor 
(, UjUj ),  the  turbulent  kinetic  energy  k=  and  the 

mean  dissipation  rate  e.  All  of  these  quantities  are  uniform  in 
space  and  evolve  in  time. 

An  essential  observation  from  experiments  and  DNS  is 
that,  after  an  initial  transient,  the  turbulence  tends  to  an  ap¬ 
proximately  self-similar  state.  The  normalized  Reynolds- 
stress  tensor 


_(uiuj) 
ij~  k 


(7) 


becomes  constant,  as  does  the  ratio  of  turbulence-to-shear 
timescales,  Sk/e,  and  hence  also  the  ratio  of  production  V  to 
dissipation  e.  The  turbulent  kinetic  energy  equation  then  dic¬ 
tates  that  k  and  e  increase  exponentially  with  time — as  is 
observed.  Thus  when  normalized  by  k  and  e,  quantities  per¬ 
taining  to  the  energy-containing  scales  of  the  turbulence  are 
self-similar.  Since  the  Reynolds  number  k2/(ev)  increases 
with  time,  small-scale  quantities  are  not  self-similar  under 
this  scaling. 

The  DNS  of  Sawford  and  Yeung4'5  are  performed  from 
the  nondimensional  time  St  =  0  until  St  =  20.  The  fluid  par¬ 
ticles  are  introduced  at  St  =  4  when  the  self-similar  state  has 
been  attained.  The  values  Sk/e  =4.83  and  V/ e  =  1.54  are  de¬ 
duced  from  the  values  of  k  and  (  it ,  u2)  from  St  =  4  until  St 
=  20;  and  the  average  value  of  the  normalized  Reynolds 
stress  tensor  over  this  time  interval  is 


C  = 


0.96 

-0.32 

0 


-0.32  0 

0.43  0 

0  0.61 


We  introduce  the  normalized  time 


(8) 


(9) 

and  the  scaled  fluctuating  velocity  following  a  fluid  particle 
A  u  (t) 

6(0-^172.  (10) 

Consistent  with  the  self-similar  state  of  the  turbulence,  we 
assume  that  u(f)  is  a  statistically  stationary  process. 


FIG.  1.  Autocorrelation  functions  pjj(s),  Eq.  (14),  from  the  DNS  data  of 
Sawford  and  Yeung  (Ref.  5). 


The  autocovariance  of  u(f)  is 

Rjj(s)  =  (ui(i)uj(t  +  s)),  (11) 

which  (in  view  of  the  assumed  stationarity)  is  independent  of 
f;  and  the  scaled  Reynolds  stress  is 

(ui(t)uj(t))  =  Cij=Rij(  0)=^p-,  (12) 


which  is  constant.  Note  that  (unlike  Ci;)  Rjj(s)  is  not  sym¬ 
metric,  although  it  has  the  property 

Rij(s)=Rji(-s).  (13) 

It  is  conventional  to  define  autocorrelation  functions  by 


Pij(s)  = 


Rn(s) 


1/2 


(14) 


(where  bracketed  suffixes  are  excluded  from  the  summation 
convention)  so  that  the  diagonal  components  of  p,;(  0)  are 
unity.  These  autocorrelation  functions  obtained  from  the 
DNS  are  shown  in  Fig.  1.  (Note  that,  by  symmetry,  p2 3 
=  P32=0-) 

The  analysis  below  shows  that  a  preferable  definition  of 
the  autocorrelations  is 


Rijis^C^R.jis),  (15) 

where  C jk  1  denotes  the  i-k  component  of  the  inverse  of  C. 
Unlike  p,; ,  R is  a  tensor,  and  at  the  origin  it  is 

Rij(0)  =  Sij.  (16) 

These  autocorrelation  functions  obtained  from  the  DNS  are 
shown  in  Fig.  2.  [There  is  a  small  inconsistency  in  the  ex¬ 
traction  of  numerical  values  from  the  DNS:  C,;  is  obtained 
as  an  average  from  St  =  4  to  St  =  20,  whereas  R  ,j(0)  is  ob¬ 
tained  at  St  =  4.  As  a  consequence,  as  may  be  seen  in  Fig.  2, 
the  numerical  values  do  not  satisfy  Eq.  (16)  exactly.] 

Based  on  R,y(s),  we  define  the  (normalized)  integral 
timescales  by 
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FIG.  2.  Autocorrelation  functions  Rjj(s),  Eq.  (15),  from  the  DNS  data  of 
Sawford  and  Yeung  (Ref.  5). 


Tij=J0  Rij(s)ds- 


The  values  deduced  from  the  DNS  data  are 


0.44 


T  = 


-0.11 


0 


-0.06  0 
0.22  0 
0  0.24 


(17) 


(18) 


By  post-multiplying  both  sides  of  this  equation  by  C_1,  we 
find  that  R(.v)  [defined  by  Eq.  (15)]  satisfies  the  same  equa¬ 
tion 


for  ,v  3=  0, 


(22) 


with  the  simple  initial  condition  R7(0)  =  I.  The  solution  to 
this  equation  (satisfying  the  initial  condition)  is15 

°°  ( —  1 )" 

R7'(s)  =  exp(-A.0=S  - ; — A"s",  for  s^0, 

n= o  n\ 

(23) 

as  may  be  verified  by  differentiating  with  respect  to  s.  It  has 
been  assumed  that  the  eigenvalues  of  A  have  positive  real 
parts,  which  is  a  sufficient  condition  for  exp(—  As)  to  con¬ 
verge  to  zero  as  s  tends  to  infinity. 

In  the  case  that  A  has  linearly  independent  eigenvectors 
the  solution  can  be  written 


Rt(s)  =  Vexp(  —  As)V  *,  for  s^0,  (24) 

and  similarly  for  R 

R(s)r=  Vexp(  — As)V_1C,  for  s^0.  (25) 

Thus  each  component  of  the  autocovariance  is  a  linear  com¬ 
bination  of  three  decaying  exponentials — decaying  because 
the  eigenvalues  are  required  to  be  positive. 

For  the  autocorrelation  timescales,  T  [Eq.  (17)]  we  ob¬ 
tain 


III.  STOCHASTIC  MODEL 

The  stochastic  model  considered  is  Eq.  (1)  written  for 
u(  /).  It  is  convenient  to  use  matrix  notation,  and  so  the  equa¬ 
tion  is  written 

du=  —Audi+BdW,  (19) 


r  cc  r  oo 

Tr=  R  T{s)ds=  exp(  —  As)ds  =  A-1.  (26) 

Jo  Jo 

The  conclusion  from  this  development  is  that  the  matrix 
of  autocorrelation  timescales  T  of  the  process  u(  t)  generated 
by  the  stochastic  model  Eq.  (19)  is  uniquely  determined  by 
the  drift  matrix  A  as 


where  (dW  dWT)  =  Idt,  with  I  being  the  identity,  and  T 
denoting  the  transpose. 

The  drift  matrix  A  is  constant  and  it  is  required  that  its 
eigenvalues  have  positive  real  parts.  The  value  of  A  deduced 
from  the  DNS  (below)  has  the  simplest  structure — real  posi¬ 
tive  eigenvalues  and  independent  eigenvectors.  In  this  case  A 
can  be  decomposed  as 

A=VAV_1,  (20) 

where  the  columns  of  V  are  the  eigenvectors  of  A,  and  A  is 
the  diagonal  matrix  of  eigenvalues. 

The  diffusion  coefficient  matrix  B  is  also  constant  and, 
without  loss  of  generality,7  we  take  it  to  be  symmetric  (B 
=  Br). 


A.  Autocorrelation  function 

It  is  readily  deduced  from  Eq.  (19)  that  the  autocovari¬ 
ance  matrix  R(,v)  [Eq.  (11)]  satisfies  the  ordinary  differential 
equation 

dRr  . 

— —  =  - ARr,  for  s^0.  (21) 

as 


T=  (A_1)r.  (27) 

This  conclusion  depends  on  the  eigenvalues  of  A  having 
positive  real  parts. 

B.  Covariance 

It  follows  from  Eq.  (19)  that  the  covariance  C  =  (uur) 
evolves  by 

dC  XT 

—  =  -AC-CAr+BBr.  (28) 

dt 

Given  that  B  is  symmetric  and  that  the  process  is  stationary, 
this  leads  to  the  relation 

B2  =  AC+CAr.  (29) 


C.  Specification  of  stochastic  model  coefficients 

Can  the  model  coefficients  A  and  B  be  chosen  so  that  the 
autocovariance  R(s)  from  the  model  matches  that  obtained 
from  DNS  of  homogeneous  turbulent  shear  flow?  Clearly  the 
answer  is  “no,”  since  the  empirical  autocovariances  will  not 
be  of  the  simple  form  implied  by  the  model — i.e.,  sums  of 
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FIG.  3.  Comparison  of  autocorrelation  functions  /),;(. S' j ,  Eq.  (14).  from  the 
DNS  data  (symbols)  and  from  the  stochastic  model  (lines)  with  coefficients 
determined  from  the  data  [Eq.  (32)  and  Eq.  (33)].  (p2 2  circles  and  solid  line; 
p33  squares  and  dashed  line.) 


three  exponentials.  Nevertheless,  the  preceding  analysis 
shows  that  A  and  B  can  be  chosen  to  match  the  covariance  C 
and  the  timescales  T.  Specifically,  given  T,  A  is  determined 
by 

A=(T~1)7'  (30) 

[see  Eq.  (27)];  then  B  is  determined  as  the  symmetric  square 
root  of 


B2  =  AC+CAr 


(31) 


[see  Eq.  (29)].  Evidently  this  specification  requires  that  T  be 
nonsingular.  An  additional  requirement  is  that  T  and  C  be 
such  that  B2  given  by  Eq.  (31)  is  positive  semi-definite. 

For  the  values  of  C  and  T  obtained  from  the  DNS  of 
homogeneous  turbulent  shear  flow,  the  values  of  A  and  B2 
obtained  from  Eq.  (30)  and  Eq.  (31)  are 


A= 


2.45  1.24  0 

0.65  4.90  0 

0  0  4.22 


and 


3.90 


B2  = 


-1.18 


0 


-1.18  0 
3.84  0 

0  5.14 


(32) 


(33) 


D.  Comparison  of  autocorrelation  functions 

Figure  3  shows  the  comparison  between  the  autocorrela¬ 
tion  functions  pjj(s)  obtained  from  DNS  compared  to  those 
from  the  model  [with  coefficients  given  by  Eq.  (32)  and  Eq. 
(33)].  Inevitably  there  are  qualitative  differences  at  the  ori¬ 
gin.  For  pu  ,  for  example,  the  DNS  value  departs  from  unity 
at  the  origin  as  1— pn(s)~s2,  whereas  the  model  departs  as 
1  — Pn(s)~|s|.  This  leads  to  the  model  values  of  pu(s) 
being  below  the  DNS  values  at  small  times;  and  then,  from 
the  matching  of  the  integral  timescales,  it  is  not  surprising 


that  at  some  later  times  the  model  value  exceeds  the  DNS 
value.  Given  these  inevitable  differences,  the  agreement  be¬ 
tween  the  model  and  the  DNS  is  as  good  as  could  be  ex¬ 
pected.  In  particular  the  model  captures  the  difference  be¬ 
tween  pn  and  the  other  two  diagonal  components  (which  are 
nearly  equal);  and  the  differences  between  p12  and  p21 . 


IV.  GENERALIZED  LANGEVIN  MODEL 


In  PDF  methods,  the  stochastic  Lagrangian  model  for 
velocity  that  is  employed  is  the  generalized  Langevin  model 
(GLM).1112'7  Applied  to  homogeneous  turbulence,  the  model 
for  u( t )  is 

d(U,) 

du  = - - - Uj  dt  +  GjjUj  dt  +  (CQe)lL  dWi ,  (34) 

"xj 

where  the  constant  C0  is  generally  ascribed  the  value  2.1. 
The  coefficient  G(J  can  depend  on  (ujUj),  e  and  d(Uj)/dXj  : 
two  particular  specifications  of  G,y  are  considered  below. 

The  transformation  of  Eq.  (34)  to  an  SDE  for  u(  1)  re¬ 
sults  in  the  general  stochastic  model,  Eq.  (19),  with  coeffi¬ 
cients 


1  IV 


A„  =  r - 1  &,■+- 


k  d(  U j)  k 


dx ; 


-  -G, 


and 


(35) 


Bu=C"2Sij. 


(36) 


Equation  (35)  can  be  rearranged  to  yield  the  value  of 
( k/s)Gjj  implied  by  the  DNS: 


k 

-G= 

e 


-2.18 

-0.65 

0 


3.59 

-4.63 

0 


0 

0 

-3.95 


(37) 


Since  B  is  found  to  be  anisotropic — as  discussed  further  in 
the  next  subsection — no  choice  of  Co  in  Eq.  (36)  yields  the 
correct  diffusion  coefficient.  Nevertheless,  the  magnitude  of 
the  diffusion  is  characterized  by 

C0=  j  trace(B2),  (38) 

the  value  of  which  deduced  from  the  DNS  is  C0  =  4.3.  By 
comparison,  the  standard  model  Eq.  (36)  yields  Co=C0 
=  2.1. 


A.  Anisotropy  of  the  diffusion  coefficient 

The  GLM,  and  also  dispersion  models,  take  the  diffusion 
coefficient  B  to  be  isotropic,  Eq.  (36).  The  reason  generally 
advanced  for  this  specification  is  consistency  with  the  Kol¬ 
mogorov  hypotheses.  For  (dimensional)  time  intervals  s  in 
the  inertial  subrange,  r^s^ikl e  (where  Tr/  is  the  Kolmog¬ 
orov  timescale),  the  Kolmogorov  hypotheses  predict  that  the 
second-order  Lagrangian  structure  function  is  isotropic  and 
linear  in  s,  i.e., 

([ui(t  +  s)-ui(t)][uj(t  +  s)-Uj(t)])  =  C0ssSij,  (39) 

where  C0  is  a  Kolmogorov  constant.  The  GLM  yields  pre¬ 
cisely  this  results  if  C0  is  taken  to  be  C0 . 
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FIG.  4.  The  scaled  tensors  B2  and  C  shows  as  ellipses  in  the  xI—x2  plane. 
The  dot-dashed  line  is  x,o CjCf/~  =  (Cii)^2,  and  the  solid  line  is  the  corre¬ 
sponding  ellipse  for  B2.  Shown  for  reference  are  dashed  lines  at  0°,  45°, 
90°,  and  135°. 


However,  the  value  of  B2  deduced  from  the  DNS  is  de¬ 
cidedly  anisotropic:  the  eigenvalues  of  B2  (which  are  all 
equal  to  C0  =  2.1  in  the  GLM)  are  found  to  be  2.69,  5.06,  and 
5.14.  It  is  possible  that  this  anisotropy  is  a  Reynolds-number 
effect,  which  vanishes  at  sufficiently  high  Reynolds  number. 
This  possibility  could  be  investigated  through  DNS  at  differ¬ 
ent  Reynolds  numbers. 

It  is  also  possible  that  the  anisotropy  in  the  deduced 
value  of  B2  persists  at  high  Reynolds  numbers,  not  because 
the  Kolmogorov  hypothesis  [Eq.  (39)]  is  incorrect,  but  be¬ 
cause  the  stochastic  Lagrangian  model,  Eq.  (19),  is  too 
simple  to  represent  the  multi-timescale  aspects  of  anisotropic 
turbulence. 

Given  the  observation  that  B2  is  anisotropic,  it  is  natural 
to  consider  modifications  to  the  GLM  to  incorporate  such 
anisotropy.  The  natural  way  to  introduce  anisotropy  in  the 
model  is  to  make  the  diffusion  coefficient  dependent  on  the 
normalized  Reynolds  stresses  C.  But  any  such  model  implies 
that  the  principal  axes  of  B2  and  C  are  aligned,  which  is  not 
supported  by  the  data.  Figure  4  shows  the  ellipses  in  the  x1 
—x2  plane  corresponding  to  the  tensors  B2  and  C.  The  mis¬ 
alignment  of  the  principal  axes  is  evident.  In  fact,  to  within 
1°,  the  minor  axis  of  B2  is  aligned  with  the  major  axis  of  the 
mean  rate-of-stain  tensor  S  (i.e.,  the  45°  line  x2=xl).  Hence 
an  anisotropic  model  for  B2  could  be  constructed  based  on  S 
that  is  consistent  with  the  DNS  data.  More  data — from  dif¬ 
ferent  flows  and  at  different  Reynolds  numbers — are  needed 
before  an  anisotropic  model  for  B2  of  any  generality  can  be 
constructed. 

B.  Simplified  Langevin  model 

In  this  and  the  next  subsection  we  examine  two  specific 
forms  of  the  generalized  Langevin  model,  corresponding  to 
particular  specifications  of  G,; . 

For  the  simplified  Langevin  model  (SLM)  considered 
here,  the  specification  is 


TABLE  I.  Values  of  the  mean  timescale  '/’  :  ]  trace(T),  the  normalized 
Reynolds  stress  CtJ ,  and  the  turbulence-to-shear  timescale  ratio  Sk/e  from 
the  DNS  of  Sawford  and  Yeung  (Ref.  5)  and  from  SLM  and  LIPM  for 
different  values  of  C0  and  a2. 


DNS 

SLM 

SLM 

LIPM 

LIPM 

C0 

- 

2.1 

3.4 

2.1 

4.4 

a2 

- 

- 

- 

3.5 

11.9 

T 

0.30 

0.43 

0.30 

0.63 

0.30 

Cn 

0.96 

1.10 

0.98 

1.02 

1.02 

C  22 

0.43 

0.45 

0.51 

0.49 

0.49 

C33 

0.61 

0.51 

0.51 

0.49 

0.49 

Cl2 

-0.32 

-0.39 

-0.34 

-0.36 

-0.36 

Sk/e 

4.83 

4.02 

4.47 

4.28 

4.28 

GU--\2  +  4CQJlSi^ 

so  that  the  matrix  A  [Eq.  (35)]  is 
X  a  0 
A=  0  X  0 
0  0  X 


with 


1  V  3 

X=2  e+4Co’ 


and 


(40) 


(41) 


(42) 


a=Sk/e. 


Evidently  all  three  eigenvalues  of  A  are  equal  to  X,  and  the 
eigenvectors  are  not  independent — two  are  equal  to  [  1  0  0]7. 


Consequently,  the  autocovariance  R(.v)  is  not  given  by  Eq. 
(25),  but  instead  the  solution  to  Eq.  (21)  is 


R(s)  =  e~ 


C  \  i  o-iC12 
^21“  °"iC22 
0 


C12 

C22 

0 


0 
0 

C33-I 


(43) 


Given  a  specified  value  of  C0  and  the  DNS  value  of 
Vie,  Eq.  (31)  can  be  solved  to  determine  the  normalized 
Reynolds  stresses  C  given  by  SLM  in  homogeneous  turbu¬ 
lent  shear  flow,  and  then  the  autocovariances  can  be  evalu¬ 
ated  from  Eq.  (43).  We  consider  two  values  of  Co:  the  stan¬ 
dard  value  Co  =  2.1;  and  the  value  Co  =  3.4  for  which  the 
SLM  timescale  X  1  matches  the  average  DNS  timescale  T 
=  j  trace  (T).  The  values  of  C  obtained  are  shown  in  Table  I. 

The  autocorrelations  p,j(s)  obtained  with  C0=3.4  are 
compared  to  the  DNS  data  in  Fig.  5.  As  expected,  the  agree¬ 
ment  is  much  better  with  the  timescales  matched  (C0  =  3.4) 
than  otherwise  (C0  =  2.1,  not  shown).  The  model  correctly 
predicts  the  equality  of  p22  and  p33  and  their  distinction  from 
Pu ,  but  the  quantitative  agreement  is  noticeably  poorer  than 
in  Fig.  3. 

The  model  predicts  a  more  substantial  difference  be¬ 
tween  p  12(s)  and  p2|(.v)  than  is  evident  in  the  DNS — a  be¬ 
havior  which  is  easily  understood.  The  only  off-diagonal 
term  [cr  in  Eq.  (41)]  enters  the  SDE  for  velocity  as 
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FIG.  5.  Comparison  of  autocorrelation  functions  p,j(s),  Eq.  (14),  from  DNS 
(symbols)  and  from  SLM  with  C0  =  3.4.  (For  SLM,  p2i  = P33  •) 


FIG.  6.  Comparison  of  autocorrelation  functions  p,;(i),  Eq.  (14),  from  DNS 
(symbols)  and  from  LIPM  with  C0  =  4.4  and  a2=  11.9  (lines).  Dashed  lines, 
P[ 2  and  p33 ;  solid  lines,  p2i  and  p22 . 


Sk 

du\= - u2dt _  (44) 

Thus  large  positive  (or  negative)  values  of  u2  tend  to  lead  to 
large  negative  (or  positive)  values  of  ux  after  a  time  lag. 
Thus  the  peak  correlation  \(u2(t)ul(t+s))\ — or  equiva¬ 
lently  the  minimum  of  p21(s) — occurs  for  a  positive  value  of 
s.  The  model  overestimates  this  effect,  because  it  takes  no 
account  of  rapid  pressure  fluctuations  which  tend  to  counter¬ 
act  the  effects  of  mean  shear. 


C.  Lagrangian  isotropization  of  production  model 
(LIPM) 

The  LIPM14'7  corresponds  closely  to  the  Launder,  Reece, 
and  Rodi16  Reynolds-stress  model.  Using  standard  values  for 
the  model  constants  /3j „ 3  and  7,  _6 ,  the  LIPM  equation  for 
G  is 


k 

e 


G=a1I+a2(b-3b2) 


5  +  5b" 


0 


(45) 


DNS  value;  see  Table  I.  As  a  consequence,  the  model  (with 
Co  =  2.1)  produces  autocorrelations  p,y(s)  (not  shown)  in 
very  poor  agreement  with  the  DNS  data. 

To  provide  a  more  meaningful  comparison,  the  constants 
C0  and  a 2  are  adjusted  to  match  the  average  timescale,  while 
leaving  the  normalized  Reynolds  stresses  the  same.  The  au¬ 
tocorrelations  given  by  LIPM  with  these  values  (C0  =  4.4, 
a2  =  1 1  -9)  are  compared  to  the  DNS  data  in  Fig.  6.  The 
agreement  is  quite  poor.  Except  at  small  times,  p22(s)  is 
incorrectly  predicted  to  be  larger  than  p33(y);  and  evidently 
the  effect  of  the  rapid  pressure  is  overpredicted  as  there  is 
little  difference  between  p12(s)  and  p2|(.v). 

This  last  point  can  be  seen  directly  in  the  matrix  A, 
which  for  LIPM  is 


A= 


2.86  1.99  0 

2.21  6.11  0 

0  0  4.48 


(48) 


The  direct  effect  of  shear  appears  in  the  1-2  component,  and 
in  SLM  the  2-1  component  is  zero.  In  A  deduced  from  the 
DNS  data  [Eq.  (32)],  A21  as  about  half  of  A 12 ;  but  for  LIPM 
A21  exceeds  A 12 . 


V.  CONCLUSIONS 


0  0  0. 
where  b  is  the  anisotropy  tensor 

b=^C-jI,  (46) 

the  standard  value  of  the  constant  a2  is  a2=  3.5,  and  the 
coefficient  a1  is  given  by 

1 1  3  \  3  V 

ai  =  -|-+-C0|  +  Y^-  +  3  a2trace(b3).  (47) 

With  the  standard  value  Co  =  2.1,  the  model  yields  rea¬ 
sonable  values  of  the  normalized  Reynolds  stresses,  but  the 
average  time  scale  7’=  j  trace! T)  is  more  than  twice  the 


As  previously  observed  by  Sawford  and  Yeung4,5  in  the 
context  of  turbulent  dispersion,  Lagrangian  data  from  DNS 
of  homogeneous  turbulence  is  valuable  in  the  development 
and  testing  of  stochastic  Lagrangian  models.  After  an  initial 
transient,  homogeneous  turbulent  shear  becomes  (approxi¬ 
mately)  self-similar,  so  that  the  appropriately  scaled  La¬ 
grangian  velocity  fluctuation  u( t)  becomes  a  statistically  sta¬ 
tionary  random  process. 

The  stochastic  Lagrangian  model  considered  for  u(  /)  is 
the  diffusion  process  Eq.  (19)  in  which  the  drift  coefficient 
depends  linearly  on  u(f)  through  the  drift  matrix  A,  and  the 
(anisotropic)  diffusion  coefficient  B  is  constant.  An  analysis 
of  this  model  shows  that  there  is  a  unique  specification  of  A 
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and  B  [Eq.  (30)  and  Eq.  (31)]  such  that  the  covariance  matrix 
C  and  the  timescale  matrix  T  match  those  obtained  from 
DNS.  The  autocorrelation  functions  predicted  by  the  model 
are  in  good  agreement  with  the  DNS  data  (except  at  small 
times).  The  model  for  u(  t)  is  a  continuous,  Gaussian,  Mar¬ 
kov  process;  and  it  is  a  significant  conclusion  that  such  a 
simple  process  provides  a  good  model  for  the  Lagrangian 
velocity  in  homogeneous  turbulent  shear  flow.  (It  is  known 
that  the  one -point  one-time  joint  PDF  of  velocity  is  jointly 
normal1  in  this  flow.) 

Contrary  to  conventional  modelling  assumptions,  if  is 
found  that  the  diffusion  coefficient  B  is  significantly  aniso¬ 
tropic.  Whether  or  not  this  is  a  low  Reynolds -number  effect 
is  an  important  question  which  can  be  addressed  in  future 
DNS  studies. 

The  magnitude  of  the  diffusion  coefficient  can  be  char¬ 
acterized  by  C0=  j  trace(B2)  and  the  value  deduced  from 
the  DNS  data  is  Co  =  4.3.  This  is  substantially  larger  than  the 
corresponding  value  C0  =  2.1  normally  used  in  PDF  models. 

There  is  evidence  that  the  appropriate  value  of  C0  de¬ 
pends  on  Reynolds  number.13'10  In  the  DNS,  the  Taylor-scale 
Reynolds  number  based  on  xi  -direction  statistics  increases 
from  Rxas40toR^s»110  during  the  course  of  the  simulation. 
Sawford  and  Yeung5  provide  an  empirical  expression  for  C0 
as  a  function  of  Rx,  which  increase  from  C0  =  3.7  at  R K 
=  40  to  C0  =  5.4  at  Rx=H0.  The  value  C0  =  4.3  deduced 
from  the  DNS  lies  within  this  range,  but  Reynolds-number 
effects  are  not  addressed  here. 

The  autocorrelation  functions  predicted  by  two  general¬ 
ized  Langevin  models  are  compared  to  the  DNS  data  in  Figs. 
5  and  6.  In  the  simplified  Fangevin  model  (SLM),  no  ac¬ 
count  is  taken  of  the  rapid  pressure  fluctuations,  and  as  a 
consequence  the  difference  between  p12 (s)  and  p2|(.v)  is 
overpredicted.  The  Lagrangian  IP  model  (LIPM) — which  in¬ 
cludes  a  model  for  the  rapid  pressure — yields  autocorrela¬ 
tions  in  poor  agreement  with  the  DNS  data. 

In  the  specification  of  both  the  drift  and  diffusion  coef¬ 
ficients,  there  is  clearly  scope  for  considerable  improvement 
in  generalized  Langevin  models. 
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